import numpy as np
from sklearn.preprocessing import PolynomialFeatures

x = [[1,2],
     [2,3]]
x = np.array(x)

pf = PolynomialFeatures()
pf.fit(x)
x = pf.transform(x)
print(x)